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Abstract. 

We have identified 335 galaxy cluster and group candidates spanning < 
z < 2, using a 4.5/mi selected sample of galaxies in a 7.25 deg 2 region in the 
Spitzer/IRAC Shallow Survey. Using full redshift probability distributions for all 
galaxies, clusters were identified as 3-dimensional overdensities using a wavelet 
algorithm. To date 12 clusters at z > 1, and over 60 at z < 0.5 have been 
spectroscopically confirmed. The mean I- [3.6] color for cluster galaxies up to 
z ~ 1 is well matched by a Zf — 3 passively evolving model. At z > 1, a 
wider range of formation histories is needed, but higher formation redshifts (i.e. 
Zf > 4 — 5) are favored for most clusters. The cluster autocorrelation function, 
measured for the first time out to z — 1.5, is found not to have evolved over 
the last 10 Gyr, in agreement with the prediction from ACDM. The average 
mass of the IRAC Shallow Cluster Survey sample, inferred from its clustering, 
is ~ 10 14 M CT . 



Introduction 



The bulk of the stellar mas s in the universe is created at 1 < z < 3 ([Dickinson et al. 
20031 : iRudnick et alJl200fih in very massive lum i nous and ultralumino us galax- 



ies (LIRGs and ULIRGs, iLe Floc'h et al.1 120051 : iPapovich et al. 
studies show that the SFR at z = 1 increases with increasing local density 



2006). Recent 



(lElbaz et all 120071: Cooper et al J 1200811. a reversal of the local SFR-density re- 
lation ( Lewis et all 20021 : Brinchmann et al. 2004). To obtain a complete cen- 
sus of massive galaxy formation at z > 1 we need to study it across a range 
of environments, not least because most massive galaxies are found in very 
rich environments. Here we describe the IRAC Shallow Cluster Survey (ISCS; 
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Eisenhardt et al. 20081 . hereafter E08), the largest, best-studied sample of stellar 



mass selected clusters at z > 1. 



2. The IRAC Shallow Cluster Survey 



The ISCS is a sample of 3 35 galaxy clusters span ning 0.1 < z < 2 in the 
IRAC Shallow Survey (ISS, lEisenhardt et~aD 120041) area of the NOAO Deep- 
Wide Field Survey (NDWFSH annuzi and DevlTl999: Jannuzi et al. in prep) in 
Bootes. The clusters are selected using a wavelet detection algorithm which 
identifies peaks in cluster probability density maps constructed from accurate 
photometric redshift probability function s for 175,431 galaxie s brighter than 



13.3/xJy at 4.5/um in a 7.25 deg 2 region (jBrodwin et al 
E08). The redshift distribution is shown in Fig. [TJ 



2006, hereafter B06, 
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Figure 1. Observed redshift distribution of 335 galaxy clusters in the ISCS. 
The curve, a fit to the distribution, is used in the clustering analysis in §4. 

The AGES survey (Kochanek et al. in prep) in Bootes provides spectroscopic 
confirmation for over 60 clusters at z < 0.5 (E08). At higher redshift, a multi- 
year Keck spectroscopic campaign has to date confirmed 12 z > 1 clusters, many 
with ~ 10 — 20 m embers, at redshifts from 1.06 to 1.41 ( Stanford et al.ll2005l : 
Elston et al.ll2006l : B06; E08). 



3. Formation of Massive Galaxies 

The extensive redshift baseline of the ISCS sample permits a test of various plau- 
sible star formation histories (E08). Fig. [2] (left) plots the average k-corrected 
I- [3.6] color for galaxies within 1 Mpc of each cluster center whose integrated 
redshift probability distribution in the range z c \ ± 0.06(1 + z c \) exceeds 0.3 (cir- 
cles; filled symbols indicate z > 1 spectroscopically confirmed clusters. The red 
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curve is not a fit, but rather a Bruzual and Chariot ( 20031 ) population synthe- 
sis model in which the stars formed in a 0.1 Gyr burst at zt = 3 and evolved 
passively thereafter. This scenario clearly describes the data better than a no- 
evolution model. Also plotted are the average k-corrected colors for galaxies 
more massive than a passively evolving L* galaxy (red boxes). These are sys- 
tematically redder than the full cluster population, showing the persistence of 
the color-magnitude relation in clusters out to z w 1.5, and suggesting the early 
existence of a mass-metallicity relation in clusters and/or providing evidence of 
downsizing. 




Figure 2. (left): Mean, k-corrcctcd I- [3.6] color vs. redshift for clusters in 
the ISCS (circles; filled symbols are confirmed z > 1 clusters). The red curve 
is a passively evolving Zf = 3 model. Red boxes are the mean colors of 
massive (L > L*) galaxies only and are systematically redder than the overall 
cluster population, right: Mean I- [3.6] color relative to the Zf = 3 PE model. 
Clusters at z ^ 1.2 formed their stars at very high redshift (zj J> 4). 

In Fig. [2] (right) the average k-corrected color of L > L* cluster galaxies is 
plotted relative to the Zf = 3 passive evolution model (horizontal line). Over- 
plotted are models with formation redshifts of zj = [1.5,2,4,5,30], as well as a 
no-evolution model. The offset between all models and the data at low redshift 
is primarily due to the redder colors of these massive galaxies discussed above. 
Nevertheless we can draw two robust conclusions from this figure regarding the 
star formation histories of z > 1 clusters. First, the variation in the average clus- 
ter colors increases markedly at z > 1. Second, we find that the stars in most 
high-redshift clusters, particularly those above z > 1.2, formed at z ~ 4 — 5 or 
even earlier. 



4. Clustering of Galaxy Clusters to 



1.5 



The large, uniformly selected ISCS cluster sample has permitted the first mea- 
surement of the galaxy cluster autocorrelation function in the first half of the 



Universe ( Brodwin et al. 20071 ). Specifically, the autocorrelation function was 
measured in two broad bins, at 0.25 < z < 0.75 and 0.75 < z < 1.5, with result- 
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ing clustering amplitudes of r = 17.40_3^q and r = lQ.MI^gg h -1 Mpc, cor- 
responding to average halo masses of logpWboo/-^©] ~ 13.91q2 an d ~ 13.81q3, 
respectively. Errors are computed via bootstrap resampling of the cluster catalog 
and are therefore rather conservative. 

These correlation lengths are quite similar to the v alues for local cluster s 
(e.g. Abell, APM, SDSS) as well as those at z « 0.5 (|Gonzalez et al.1 120021 ) . 



This constancy with redshift is a key prediction of ACDM, as shown in Fig. 
[3] (left). Measurements of the cluster clustering amplitude, ro, are plotted vs. 
intercluster distance, d c , for several different surveys spanning redshifts up to 
z m 0.5. Overplotted are the present measurements, at effective redshifts of 0.5 
and 1.0 (filled symbols). The h ashed regions shows the prediction from a high- 
resolution numerical simulation (jYounger et al. 20051 ) that over < z < 1.5 very 
little evolution is expected, in agreement with our results. 
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Figure 3. Left: r vs. d c for the present sample (filled circles) along with 
several z ^ 0.5 measurements from the literature (jBahcall et al.l [2003L and 
reference s therein). The hash ed region shows the ACDM prediction over < 
z < 1.5 ([Younger et al.ll2005h . Right: Comoving correlation lengths for the 
ISCS clusters, other cluster samples at low-z, and high-z, highly clustered 
galaxy populations. The filled region (dotted lines) illustrates the predicted 
ro evolution in a biased structure formation model for our z e g = 0.53 (0.97) 
cluster samples. The dashed lines represent simple stable clustering models. 



In Figure [3] (right) we plot a compilation of recent clustering amplitudes 
for variou s cluster surveys, as well a s for h ighly clustered galaxy populations. 
Following iMoustakas and Somervi 13 (I2002I ). we overplot the halo conserving 
model of iFrvl (11996) normalized to our two measurements in order to explore 
possible evolutionary connections with structures at other redshifts. The shaded 
area (dotted lines) shows the 1 a region for the z c s = 0.53 (0.97) measurement. In 
this model the ISCS clusters will evolve into typical present-day massive clusters, 
such as those in the SDSS, APM or Abell surveys. In the stable-clustering 



pictu re, in which clustering is fixed in physical coordinates (jGroth and Peebles! 
19771 . dashed lines), the z e s = 0.97 ISCS clusters grow into the most massive 
clusters in the local Universe, typically identified in X-ray surveys. 
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Most of the plotted high redshift galaxy clustering measurements are rather 
uncertain due to both small number statistics and poorly known redshift distri- 
butions. Several authors also find clustering amplitud es for galaxy popula tions 
which are inconsistent with their space densities (e.g. iQuadri et al.l 20071 ). In- 
deed, only obscured ULIRGs have space densities similar to the present cluster 
sample, a prerequisite for drawing evolutionary connections from these particular 
models. Their clustering is consistent with that of the ISCS clusters, providing a 
measure of support for recent studies (e.g. iMagliocchetti et alJfeOQTl ) indicating 
that they may be associated with, or progenitors of, groups or low-mass clusters. 
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